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ABSTRACT 


The  purpose  of  this  research  is  to  explore  the  need  for 
time-dependent  sea-surface  temperatures  in  atmospheric  model 
predictions  to  10  days.  Six  and  nine-layer  versions  of  the 
Navy  Operational  Global  Prediction  System  (NOGAPS)  are  used 
in  this  study.  Control  forecasts  were  made  in  which  the 
sea-surface  temperature  (5ST)  is  nxed  in  time.  Test  hind- 
casts  were  made  in  which  the  SST  was  updated  at  each  time 
step  of  the  atmospheric  model  using  interpolations  of 
Id-hourly  SST  analyses.  The  10-day  predictions  are  compared 
to  determine  any  improvement  or  degradation  due  to  the  time- 
dependent  SST.  Two  cases  are  analyzed,  one  during  Nov  m  her 
1933  and  another  during  April  1984.  Use  of  the  time- 
dependent  SSI's  resulted  in  significant  cnanges  in  the  fore¬ 
cast  fields  of  surface  heat  fluxes  and  precipitation  winch 
were  physically  consistent  with  tne  SST  trenl.  Analysis  of 
13  storm  forecasts  revealed  significant  changes  of  storm 
track,  duration  or  cyciogenesis  in  only  4  cases.  Three  of 
these  cases  were  forecast  by  the  nine-layer  version  of 
NO  GAPS  during  the  April  period  and  one  case  was  forecast  by 
the  six-layer  NOGAPS  during  the  November  period. 
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I.  IHTRODOCTIQN 


Atmosphere  and  ocean  interaction  is  recognized  as  an 
important  component  of  the  atmospheric  circulation  on  time 
scales  of  a  month  and  longer-  Air-sea  interactions  or. 
shorter  time  scales,  such  as  tnose  over  wnioh  numerical 
weather  prediction  (NWP)  models  are  integrated,  are  less 
well  understood-  Operational  atmospheric  prediction  models 
generally  employ  one-way  ocean  to  air  influence  derived  from 
a  time-independent  sea-surface  temperature  (S3!)  field.  Tne 
guaiitv  of  forecasts  from  such  models  declines  rapidly  when 
extended  beyond  five  or  six  days.  Many  factors  contribute 
to  this  observed  degradation:  inaccurate  specification  of 
the  initial  conditions;  inadequate  horizontal  and  vertical 
resolutions;  inaccurate  or  simplistic  parame te ri za ti ons  of 
atmospheric  processes;  and  lack  of  two-way  interaction 
between  the  atmosphere  and  ocean. 

Sea-surface  temperature  changes  alter  the  noundary 
forcing  of  atmospheric  prediction  models  through  the  surface 
fluxes  of  sensible  heat,  latent  heat  and  momentum.  The  heat 
fluxes  are  functions  of  the  air-sea  temperature  difference 
and  planetary  boundary  layer  (PEL)  conditions.  These  fluxes 
can  also  affect  the  radiation  balance  tnrougn  par aaeteriza- 
tions  of  cloud  cover  and  water  vapor  content.  The  magnitude 
of  the  atmospheric  response  to  changing  SST's  should  be 
related  to  the  magnitude  and  locations  of  SST  changes.  The 
nonlinear  character  of  these  relationships  requires  the  use 
of  a  global  atmospheric  prediction  model  with  sophisticated 
parameterizations  of  the  interactions  between  surface 
fluxes,  clouds,  precipitation  and  radiative  fluxes. 

This  work  is  the  second  in  a  series  of  case  studies 
designed  to  study  the  necessity  and  feasibility  of  coupling 


at  losphenc  and  oceanic  models.  The  first  case  study 
(Fanelli,  1984)  examined  the  response  of  10-dav  atmospheric 
forecasts  to  olserved  SSI  changes  introduced  tnrough  simu¬ 
lated  n on- synchronous  coupling.  Trie  temperatures  at  the  sea 
surface  were  updated  ty  replacing  the  initial  SST's  with  the 
Fleet  Numerical  Oceanography  Center  (FNOC)  3ST  analyses 
valid  at  12-hour  intervals  during  the  forecast  period.  In 
this  second  step,  simulated  syncnronous  coupling  is  used  to 
introduce  olserved  SST's  into  the  atmospheric  model.  Two 
cases,  one  each  in  the  spring  anu  autumn  transition  periods, 
are  investigated.  While  a  global  model  x s  used,  the  anal¬ 
ysis  is  restricted  to  the  northern  hemisphere  mid- lat it ude 
oceanic  areas  and  the  effects  of  changing  SST's  on  the 
atmospheric  model  10-day  predictions. 


II.  B ACKGEOJND 

Currtrt  operational  NWP  models  fall  far  short  cf  the 
approximat eiy  15-day  theoretical  limit  of  atmospheric 
predictability  (Fosmond  et  al.  1  983).  Many  methods  are 
being  proposed  and  tested  to  incrementally  improve  existing 
models.  These  include  improved  numerical  differencing 
schemes,  spectral  model  formulation,  increased  vertical  and 
horizontal  resolution,  improved  initialization  techniques, 
and  more  accurate  parameterization  of  diabatic  processes. 
While  efforts  in  these  and  other  areas  will  ultimately  be 
r«_guired  in  the  process  of  NKP  model  improvement,  the 
present  study  will  concentrate  on  the  effects  of  diabatic 
processes  in  NWP  models. 

Diabatic  processes  become  more  important  in  determining 
atmospheric  circulations  at  long  time  scales.  At  shorter 
time  scales  {1-3  days),  dynamic  processes  seem  adequate  as 
the  primary  forcing  in  NWP  models.  Diabatic  processes  begin 
to  become  important  in  medium  time  scales  (5-15  days),  and 
tend  to  dominate  the  long  time  scale  (15-30  days  and  beyond) 
forecast  problem.  Because  the  ocean  is  a  major  heat  and 
moisture  source,  the  air-sea  fluxes  are  very  important  in 
modeling  diabatic  processes. 

Latent  and  sensible  heat  fluxes  across  the  sea  surface 
are  modelled  as  functions  of  air-sea  temperature  differ¬ 
ences.  Wnile  Sandgathe  (1981)  concluded  that  maritime 
cyclogenesis  required  accurate  specification  of  air-sea 
fluxes,  no  operational  model  includes  time-dependent  SST. 
The  European  Center  for  Medium- range  Weather  Forecasts 
(ECMFF)  and  the  'Jnited  States  National  Meteorological  Center 
pNMC)  ust  tne  NNC  weekly  global  SST  analysis  in  initializa¬ 
tion.  The  ENOC  NOGAPS  model  is  initialized  by  12-hourly  SST 


analysts  produced  by  the  T her mody namic  Ocean  Prediction  and 
Expanded  Ocean  Thermal  Structure  Systems  (TCPS-EOTS) 
described  by  Clancy  and  Pollack.  (  1  983).  v:ithout  a  time- 
dependent  specification  of  SST,  the  model  representation  of 
air-sea  interface  fluxes  and  cyclogenesis  is  bound  to  dete¬ 
riorate  with  increasing  forecast  time. 

Significant  SST  changes  can  occur  during  5-15  day 
periods.  Atmospheric  forcing  causes  large  SST  changes 
uucing  the  spring  and  autumn  transition  periods  (Camp  and 
Elsberry,  1978;  Elsberry  and  Camp,  1978;  Elsberry  and  Raney, 
1978).  Curing  these  periods,  warm  shallow  oceanic  mixed 
layers  are  mechanically  mixed,  deepen  and  cool.  Winter  SST 
changes  are  generally  smaller  than  transition  changes 
because  of  a  deeper  ocean  mixed  layer.  Summer  SST  changes 
art-  smaller  due  to  decreased  mechanical  mixing  (fewer  and 
less  intense  storms).  The  intensity  and  orientation  cf  SST 
gradients  associated  with  western  boundary  currents  are 
important  in  determining  cyclogenesis  in  the  western  ocean 
lasir.s  (Sanders  and  Gyakum,  1989).  Significant  SST 
increases  in  tropical/equatorial  regions  are  manifest  in 
local  enhancements  ir.  the  convective  cloud  amounts  on  short 
(1-2  days)  time  scales. 

The  specification  of  a  time-varying  SST  to  a  NKP  model 
would  require  a  corresponding  oceanic  prediction  model.  The 
inclusion  of  atmospheric  forcing  influences  on  SST  changes 
would  require  an  input  to  the  oceanic  model  from  the  atmos¬ 
pheric  model.  The  sensitivity  of  both  models  to  the  input 
from  the  other  must  be  well  understood  before  truly  coupled 
atmospheric-oceanic  models  become  operational. 


A.  COOPLIHG  SCHEHES 

Three  methods  of  model  coupling  were  described  by 
Ei sherry  et  al.  (1982)  : 

1.  minimal  feedback; 

2.  n on- sy nchronous  coupling;  and 

2-  s y nchronous  coupling. 

Minimal  feedback  (Fig.  la)  is  tne  metnod  used  by  current 
operational  models.  It  involves  sp  ecir  ication  of  SST  at  the 
initial  time  a Dd  no  time-dependence  of  the  SST  field. 
Hon-synchronous  coupling  (Fig.  1b)  involves  SST  input  from 
an  essentially  independent  ocean  prediction  model. 
Synchronous  coupling  (Fig.  1c)  involves  concurrently  running 
oceanic  and  atmospheric  models  wnich  provide  boundary 
forcing  to  each  other. 

An  understanding  of  the  effects  of  more  sophisticated 
coupling  schemes  on  NrfP  models  is  required  before  these 
schemes  are  implemented  in  operational  models.  The  addi¬ 
tional  degrees  of  freedom  afforded  by  the  new  time  dependen¬ 
cies  in  boundary  forcing  can  possibly  degrade  model  forecast 
skill,  especially  that  achieved  tnrough  "tuning"  of  parame¬ 
terized  processes.  Atmospheric  model  biases  may  result  from 
biases  in  tbe  oceanic  prediction  models  or  improper  coupling 
between  models.  A  large  number  of  model  verification 
studies  will  be  required  to  identify  these  biases  with  any 
statistically  significant  degree  of  certainty. 

B.  PREVIOUS  STUDIES 

Many  previous  studies  have  concentrated  on  long  time 
scale  (cli mat  a  1 ogi ca 1)  effects  of  SST  changes  on  the  atmos¬ 
phere.  A  review  of  these  studies  is  contained  in  Elsterry 
e  t  a  1 .  (18  82). 
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Arpe  (1931)  examined  the  effect  of  artificial  S5T  anoma¬ 
lies  on  the  performance  of  the  EC MW F  model.  He  shoved  model 
s>.-r.si t i vit y  at  forecast  times  of  six  days  and  increased 
cyci  oaer.esis  with  large-scale  positive  351  anomalies. 
Fanelli  (1984)  examined  the  effects  of  observed  SST  devia¬ 
tions  (with  respect  to  the  initial  conditions)  on  the  NOGAPS 
codui.  clear  differences  in  model  nehavior  were  noticed 
over  a  10-day  forecast.  The  model  employing  the  time- 
impendent  SST  seemed  to  show  improved  skill  in  predicting 
mil-latitude  maritime  storm 


lifetimes. 


III.  PROCEDURE 


The  experimental  procedure  was  similar  to  that  of 
Fanelli  (1984) .  Two  versions  or  NOCAPS  were  used  in  tnis 
research.  In  each  of  two  cases,  two  N03APS  integrations  to 
10  days  forecast  time  were  made.  The  first  of  the  integra¬ 
tions  was  the  control  run  in  which  the  SST  was  held 
constant.  In  the  second  integration,  designated  the  SST 
run,  global  SSI's  were  updated  at  every  time  step  (4 
minutes)  of  the  model.  The  model  atmospheric  responses  of 
each  forecast  were  analyzed  for  differences.  The  analysis 
was  restricted  to  the  following  fields: 

a.  surface  pressure; 

b.  sea-surface  temperature; 

c.  surface  sensible  heat  flux; 

d.  surface  latent  neat  flux; 

e.  surface  total  heat  flux  (includes  radiation); 

f.  convective  precipitation;  and 

g.  large-scale  precipitation. 

The  NOGAPS  model  was  made  available  by  Dr.  T.  Rosmond  of 
the  Naval  Environmental  Prediction  Research  Facility 
(NEPRF)  .  The  general  caa ra cte  r istics  of  the  model  were 
described  by  Ranelli  (1984).  During  the  course  of  this 
experiment,  major  changes  were  made  to  NOGAPS.  The  November 
1983  case  study  was  run  with  a  six-layer  version  of  NOGAPS 
very  similar  to  that  used  by  Ranelli  (1984).  The  April  1984 
case  was  rur.  with  the  new  nine-layer  version  of  NOGAPS.  In 
alditior.  tc  vertical  resolution,  tne  two  model  versions 
differed  in  the  method  by  which  diabatic  heating  was 
applied.  In  the  earlier  version,  tne  diaratic  processes 
were  calculated  ar.d  applied  once  every  Id  model  time  steps 


(4  0  minutes).  The  new  version  also  calculates  the  diabatics 
every  10  model  tine  ste)  s  but  applies  tuis  heating  uniformly 
at  each  time  step.  Consequently,  the  diaratic  effects  will 
he  a  i  plied  smoothly  with  less  shock  between  the  time  steps 
ir.  which  diabatic  effects  are  and  are  not  introduced.  The 
effects  of  time- dependent  SST  will  also  be  introduced  more 
smoothly  in  the  new  NOGAPS.  Since  additional  model  levels 
are  included,  the  coupling  between  the  planetary  boundary 
layer  and  the  free  atmosphere  should  also  be  enhanced. 

Observed  sea-surface  temperatures  were  use!  to  simulate 
tna  input  from  an  ocean  prediction  model.  The  SST  analyses 
are  produced  twice  daily  for  valid  times  of  0000  GMT  and 
1200  GMT.  These  analyses  were  used  to  create  fields  of 
12-hour  forward  difference  SST  change.  NOGAPS  was  modified 
to  input  these  12-hour  forward  difference  SST  changes  and 
apply  an  incremental  portion  of  the  SST  change  at  each  time 
step.  This  process  simulates  a  synchronous  coupling  of 
NOSAPS  to  a  "perfect  prog"  ocean  model. 

The  NOGAPS  initial  conditions  used  for  the  case  studies 
were  identical  to  those  used  for  the  operational  forecasts 
with  tne  same  starting  times.  The  case  starting  times  were: 

a.  1200  GMT  31  Gctober  1963;  and 

b.  1200  GMT  21  April  1984. 

Ti.e  model  output  parameters  were  analyzed  by  three 
methods.  Forecast  and  analyzed  vortices  were  objectively 
tracked  using  the  Systematic  Errors  Identification  System 
(SETS)  developed  by  Brody,  et  al.  (1984)  at  NEPEF.  SEIS  is 
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cast  differences  and  the  time- dependent  SST. 
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IV.  HEAT  FL OX  ANALYSIS 


The  introduction  of  time-dependent  sea-surface  tempera¬ 
tures  into  the  NOGAPS  model  results  in  immediate  changes  in 
the  surface  fluxes  of  sensible  and  latent  heat.  These 
fluxes  in  KG3APS  are  parameterized  according  to  Peardorff 
(  1972).  Because  the  PEL  and  convective  parameterization 
schemes  in  NOGAPS  are  intimately  coupled,  the  changes  in 
air-sea  fluxes  may  change  the  diabatic  forcing  throughout 
the  troposphere. 

A.  SEA-SURFACE  TEMPEBATUEE  CHANGES 

In  both  the  November  and  April  cases,  significant  SST 
cnanges  evolved  luring  tne  10-day  forecast  period.  The  SSI 
difference  field  (day  10  -  initial  conditions)  for  the 
November  case,  is  shown  in  Fig.  2.  Lace-scale  negative  SST 
changes  are  evident  in  the  western  Pacific  Ocean,  Sea  of 
Okhotsk,  eastern  Pacific  Ocean,  and  mid- Atlantic  Ocean. 
Ti.ere  weie  small  areas  of  weak  positive  SSI  change  to  the 
west  of  the  3ritish  Isles  and  to  the  east  of  the  nortn- 
eastern  United  States.  The  SST  difference  fields  for  the 
ALrii  case  are  shown  in  Fig.  3.  Fhile  both  positive  and 
negative  SSI  changes  are  shown,  the  overall  mean  SST  differ¬ 
ence  is  positive.  Large  positive  differences  are  found  in 
the  Sea  of  Japan,  Yellow  Sea,  eastern  Pacific  and  western 
Allan- ic  Oceans.  Negative  differences  (decreasing  SST  with 
time)  occur  in  the  western  Pacific  (east  of  Japan)  and  north 
Atlantic  (south  of  Greenland).  The  evolution  of  the  field 
mean  SST  differences  for  eaci.  ocean  oasin  in  botn  cases  is 
shown  in  Fig.  4.  The  10-day  trends  in  toe  mean  SSI  differ¬ 
ences  are  essentially  linear.  Correlation  coefficients  of 
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linear  regressions  on  the  field-mean  SST  differences  were 
lotween  0.86'  and  0.99.  Some  of  this  trend  is  due  to  the 
actual  SST  observations.  In  addition,  the  TGPS-EGTS 
forecast/analysi s  scheme  includes  a  time-dependent  effect  in 
tne  first  guess  and  a  reversion  to  climatology  in  data- 
sparse  areas.  Each  of  these  effects  will  also  contribute  to 
a  general  trend. 

B.  SURFACE  HEAT  FLUX  CHANGES 

As  in  Fanelii  (1984),  the  surface  heat  fluxes  were 
responsive  to  the  SST  changes.  However,  the  spatial 
patterns  of  heat  flux  differences  were  not  always  related 
directly  to  the  patterns  of  SST  difference.  The  heat  flux 
difference  patterns  often  took  the  form  of  couplets  of  posi¬ 
tive  and  negative  differences. 

1 .  Sensible  Heat  Flux 

The  field  mean  heat  flux  differences  did  show  trends 
similar  to  the  SST  differences.  Tne  12-hour  fields  of 
sensible  heat  flux  differences  were  extremely  noisy. 
Ftdtures  in  the  difference  fields  were  rarely  stable  for  43 
hours.  The  largest  differences,  both  positive  and  negative, 
were  observed  in  the  western  ocean  basins  in  the  November 
case  and  in  the  Pacific  Ocean  south  of  30  N  for  the  April 
case. 

The  development  of  the  mean  sensible  heat  flux 
differences  is  shown  in  Figs.  5  tnrough  8.  While  the  field- 
c.ean  differences  reach  5-20  percent  of  the  control  forecast 
mean  sensible  heat  flux,  peak  values  in  tne  difference 
fields  exceed  50  percent  of  peak  values  in  the  control  fore¬ 
cast  fields.  Peak  differences  exceed  the  mean  control 
fluxes  by  factors  of  2-15.  Correlation  coefficients  of 
linear  regressions  of  tne  field-mean  differences  were  from 
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' .  iiiiLUL  y.eut  Flux 

The  latent  ..cat  flux  iifrerence  fields  showed  less 
1.01 : -  t  and  core  stability  that,  those  of  sensible  heat  fiux. 
A  direct  correlation  of  the  latent  flux  dirf arence  fields 
ui.d  FST  difference  fields  was  not  obvious.  Representative 
latent  heat  fiux  difference  fields  are  snown  in  Fig.  9.  In 
t ;,c  November  case,  strong  positive-negative  couplets  devel- 
opei  in  the  western  Atlantic  region.  The  Pacific  Ocean 
areas  south  of  30  N  exhibited  the  largest  differences  in  the 
April  case.  Analysis  of  the  field  statistics  (Figs.  10 
through  13)  shows  mean  differences  reaching  5-12  percent  of 
the  control  means.  Peak  values  in  the  difference  fields 
reach  25-75  percent  of  the  peak  values  of  the  corresponding 
control  latent  heat  flux  fields  and  exceed  field-mean  fluxes 
in  the  control  run  by  factors  of  2-6.  Linear  regression  of 
ti.e  field-means  produced  correlation  coefficients  of 
0.75-0.86  in  the  November  case  and  only  0.22-0.60  in  the 
April  cast.  The  differences  in  latent  heat  fluxes  after  10 
days  during  November  are  at  least  10  times  as  large  as  the 
differences  in  sensible  heat  fluxes,  and  represent  consider¬ 
able  reductions  in  evaporation  over  the  entire  field. 

3 .  Total  Heat  Flux 

NGoAPS  total  heat  flux  includes  the  sensible  and 
latent  fluxes  )  lus  the  net  radiative  fiux  at  the  surface. 
T..£  radiation  is  affected  by  model  cloudiness.  Fields  of 
t  o t  a  n eat  .lux  d i  fference  tend  to  be  correlated  with  fields 
of  latent  heat  fiux  difference.  Field  statistics  (Figs.  19 
through  17)  snowed  mean  differences  of  8-20  percent  over  the 
K-lciy  forecast  periods.  Peak  differences  at  specific  times 
were  40-  100  percent  of  the  control  peak  fluxes  and  exceeded 
control  fief.d-r.eans  by  factors  of  2-10.  linear  regression 
of  the  field-means  ~tr  educed  correlation  coef  f  icier.ts  of 
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0.44-0.  c4  ii.  the  November  case  and  0.  1u-G. 4b  in  the  April 
case.  Especially  ii.  the  November  case,  these  dirferer.ces 
atter  10  days  represent  significant  changes  in  the  surface 
heat  budget  terms  which  would  be  used  to  drive  the  ocean 
model . 


C.  SOHMARI 

The  introduction  of  time-dependent  sea-surface  tempera¬ 
tures  resulted  in  significant  field-mean  sensible,  latent 
and  total  heat  flux  differences  of  up  to  20  percent  over 
10-day  forecast  periods.  Maxima  in  the  flux  differences 
between  the  551  and  control  runs  after  10  days  were  of  the 
same  order  of  magnitude  as  pea/c  values  in  the  control  flux 
fields.  Thus,  the  surface  fluxes  at  individual  points  may 
differ  significantly  from  the  control  run  when  time- 
dependent  SSI's  are  specified.  Stable  features  in  the 
difference  fields  of  latent  and  total  heat  flux  were 
Oi served  to  intensify  and  move  consistently  during  the 
second  naif  of  the  forecast  period.  The  largest  changes 
were  noted  in  the  western  Atlantic  in  tne  November  case  and 
in  the  Pacific  (south  of  33  N)  for  the  April  case. 

Tne  forecast  differences  in  heat  fluxes  are  of  such 
magnitude  that  they  can  be  expected  to  change  the  nature  of 
the  a tmosr heric  (and  oceanic)  forecast  at  some  future  time. 
Local  changes  in  the  heat  fluxes  could  also  be  expected  to 
seriously  affect  model-derived  forecasts,  such  as  those 
produced  by  model-output  statistics  and  systems  similar  to 
tne  Navy  Operational  Local  Atmospheric  Prediction  System 
(I.'OLAPS)  ,  since  these  products  depend  on  the  predicted 
surface  fluxes. 


Ti.e  oranges  in  the  surface  neat  fluxes  were  greater  than 
those  on-served  !  y  Ranelli  (1934;.  Tnree  Key  differences  in 
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V.  PRECIPITATION  ANALYSIS 


Changes  in  the  surface  heat  fluxes  induced  by  time- 
dependent  sea-surface  temperatures  can  be  expected  to  affect 
convective  and  large-scale  precipitation  forecasts.  A  study 
by  Bosse  (1984)  showed  that  surface  heat  fluxes  could  be 
related  to  NOGAPS  large-scale  precipitation  fields  through 
layer  cloud  instability. 

A.  CUHULUS  PRECIPITATION 

The  development  of  the  cumulus  precipitation  difference 
fields  during  tne  10-day  forecasts  is  shown  in  Figs.  18 
through  21.  Over  both  the  North  Pacicic  and  North  Atlantic, 
there  are  25  to  50  percent  decreases  in  mean  cumulus  precip¬ 
itation  during  the  November  forecast.  In  the  April  fore¬ 
cast,  mean  cumulus  precipitation  over  the  Pacific  decreased 
almost  15  percent  while  cumulus  precipitation  over  the 
Atlantic  increased  20  percent.  Peak  difference  values 
exceeded  100  percent  of  peak  at  some  time  during  each  of  the 
forecasts.  The  decrease  in  cumulus  precipitation  during  the 
November  period  seems  logical  in  view  of  the  decreasing  mean 
sea-surface  temperatures  and  latent  neat  fluxes.  However, 
the  magnitude  of  the  changes  is  surprisingly  large. 
Correlation  coefficients  of  0.74  were  obtained  in  linear 
regression  of  both  the  Atlantic  and  Pacific  near,  cumulus 
precipitation  differences.  The  magnitudes  of  the  SST  and 
latent  heat  flux  field-mean  differences  were  much  less  in 
ti.e  April  case  than  in  the  November  case.  The  variance  of 
Lotn  mean  SST  change  and  mean  flux  change  was  greater  in  the 
Aj.  ril  case  than  in  the  November  case  due  to  tne  presence  of 
areas  of  strong  positive  ar.a  negative  SST  change  in  both 
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:  .ii  i:.:  .  This  cay  lain  tne  ltss  consistent  changes 
iii  cuialui  j  tix'ij  i tut  ion  curing  the  April  forecast, 
roru  lutiui.  cjvf  iicui.ts  for  the  linear  regression  of  the 
iti:.  cur.uiu;'  iifferences  luring  the  April  period  were  O.fal 
in  tne  i  acif ic  and  0.007  in  the  Atlantic. 

b.  LARGE-SCALE  PRECIPITATION 

The  development  of  the  large-scale  precipitation  differ¬ 
ence  fields  is  shown  in  Figs.  22  through  25.  Mean  differ¬ 
ences  or  large-scale  precipitation  are  less  than  10  percent 
except  for  tae  April  Atlantic  at  the  end  of  the  forecast 
period.  Peak  differences  range  from  75  to  100  percent  of 
peak  control  values. 


C.  SUMMARY 

Significant  mean  and  peak  precipitation  differences  were 
of  served  during  the  10-day  forecasts.  As  with  the  heat  flux 
differences,  these  changes  in  model  precipitation  grow  with 
time.  Long  range  model-derived  products  such  as  model- 
output  statistics  will  be  affected  by  these  differences. 


VI.  ANALYSIS  OF  STORM  TRACKS 


The  S El S  program  (Erody,  et  al.,  1984)  was  used  to  objec¬ 
tively  analyze  the  position,  size,  shape  and  amplitude  of 
low  pressure  systems  in  the  control  and  SST  forecasts  and 
tne  analyses.  The  analysis  was  performed  on  surface  pres¬ 
sure  only. 

The  SETS  storm  parameters  are  defined  in  a  manner 
similar  to  that  of  Williamson  (1981).  Amplitude  is  the 
difference  between  the  storm  central  pressure  and  the 
latit udinaily-averaged  surface  pressure.  Radius  is  the  size 
of  the  storm  expressed  in  terms  of  tne  spacing  of  the  FNOC 
Northern  Hemisphere  polar  stereographic  grid.  Ellipticity 
is  the  square  root  of  the  ratio  of  the  semi-major  and  semi¬ 
minor  axes  of  the  ellipse  which  best  fits  the  analyzed  storm 
snape.  Angle  is  the  angle  between  the  x  axis  and  the  storm 
ellipse  major  axis,  measured  counter  clockwise  from  the  x 
axis . 

In  the  SETS  program,  the  forecast  io#  pressure  systems 
die  analyzed  first,  then  low  pressure  systems  in  the  anal¬ 
yses,  then  forecast  and  analyzed  systems  are  matcned 
according  to  set  criteria.  The  nominal  5EIS  matching 
criteria  are  plus  or  minus  2.5  grid  points  in  both  the  x  and 
y  directions  and  plus  or  minus  U  hours  in  time.  5EIS  is 
routinely  used  to  evaluate  60  hoar  forecasts  and  the 
matching  criteria  are  optimized  for  this  period.  As  a 
result,  the  vortex  tracking  and  description  portions  of  SE1S 
perform  satisfactorily  throughout  the  10-day  forecasts,  but 
the  automatic  matching  portion  of  SEI5  necom.es  ineffective 
after  five  days.  This  is  primarily  a  function  of  the  overly 
restrictive  time  matching  criteria.  Inspection  of  the  SETS 
out;  ut  indicated  that  several  good  storm  traci  matches  in 


the  5-10  day  portion  of  the  forecast  period  could  be  made  if 
the  time  matching  criteria  were  increased  to  plus  or  minus 
48  hours.  Manual  matching  with  tne  standard  distance 
criteria  and  48  hour  time  criteria  was  used  in  the  5-10  day 
portion  of  the  forecasts. 

Labels  were  assigned  to  the  matcned  storm  tracks  for 
discussion  and  reference  purposes.  Tne  storm  labels  consist 
of  3  characters.  The  first  cnaracter  is  the  first  letter  of 
the  forecast  month  (N  for  November,  A  for  April)  ;  the  second 
character  is  the  first  letter  of  the  ocean  basin  (A  for 
Atlantic,  P  for  Pacific) ;  and  tne  third  character  is  a 
numeral  assigned  on  the  basis  of  chronological  order  of 
storm  identification  within  a  particular  forecast  and  ocean 
basin. 

A.  STORMS  DURING  NOVEMBER 

The  November  control  and  S3T  forecasts  were  made  with 
tne  six-level  version  of  NOGAPS.  A  total  of  15  storms  with 
lifetimes  greater  than  36  hours  were  tracked  in  the  surface 
pressure  analyses.  Fifteen  storms  were  present  in  the 
control  forecast  while  16  were  tracked  in  tne  SST  forecast. 
Matches  were  obtained  for  eight  of  the  15  analyzed  storm 
tracks. 

relatively  minor  differences  between  the  control  and  SST 
forecast  storm  positions  and  parameters  were  noted.  Tnese 
differences  grew  with  time.  The  traces  for  storm  N A 5,  which 
si. owed  the  largest  track  and  parameter  differences,  are 
shown  in  Fig.  26.  This  storm  was  first  identified  at  72 
hours  in  the  analyses  and  at  56  hours  in  noth  .forecasts. 
Tne  analysis  and  control  forecast  storm  lifetimes  of  120 
hours  are  identical  while  the  SST  storm  lifetime  was  24 
hours  longer.  The  forecast  storm  tracks  diverge  after  the 
third  forecast  hut  remain  almost  naraiiel.  The  control  run 


storm  rtiLdii.s  q  uasi-stationary  north  of  Great  Britain  four 
days  after  first  detection,  while  the  SST  run  storm  stalls 
for  only  If  hours  in  the  same  region  and  tuer.  proceeds  to 
the  northeast.  The  cyclogenesis  and  early  storm  tracks  of 
NA 5  are  in  the  region  of  maximum  SST  change,  but  the  late 
tracks  and  cyclolysis  positions  are  in  the  region  of  minimum 
SSI  change.  The  SEIS  parameters  for  storm  N A 5  are  shown  in 
Fig.  27.  The  forecast  amplitude  and  period  of  maximum 
radius  lag  the  analyzed  parameters  by  48  to  b3  hours.  The 
control  rur.  was  more  accurate  than  the  SSI  run  in  this 
particular  case,  although  the  differences  are  much  smaller 
tnar.  the  departures  of  the  forecasts  from  the  analyses.  The 
I.'A5  forecast  was  the  only  one  of  the  eight  November  storms 
for  which  one  forecast  could  be  judged  more  accurate  tnan 
the  other.  The  descriptions,  tracks  and  parameters  for  the 
other  seven  November  storms  are  contained  in  Appendix  A. 


B.  STORMS  DURING  APRIL 


The  April  control  and  SST  forecasts  were  made  with  the 
nine-layer  version  of  NOGAPS.  A  total  of  11  storms  with 
lifetimes  greater  than  or  equal  to  36  hours  were  tracked  in 
the  surface  pressure  analyses.  Fourteen  were  tracked  in  the 
control  forecasts,  while  16  were  present  in  the  SST  fore¬ 
casts.  Matches  were  obtained  for  seven  analyzed  storm 
tracks.  Significant  differences  between  control  and  SST 
forecast  tracks  and  parameters  were  found  for  tnree  of  the 
seven  matched  tracks.  Those  matched  tracks  witn  significant 
differences  are  discussed  here.  The  remainder  of  the  April 
storm  tracks  are  discussed  in  Appendix  A. 

1 .  Storm  A  A3 


The  tracks 


and  SEIS  parameters  for  storm  A A  3  are 
Tnis  storm  was  first  detected  in 


s..owr.  in  F 


15  and 


the  analyses  at  84  hours  and  was  forecast  only  by  the  SST 
run  at  108  hours.  The  forecast  storm  was  weaker  (5  mb)  and 
shorter-lived  (36  vice  168  hours)  than  the  actual  storm,  but 
the  cyclogenesis  position  and  initial  track  were  well 
predicted.  This  was  the  only  case  for  which  only  one  of  the 
forecast  models  predicted  a  matched  storm  track  when  the 
other  forecast  model  did  not.  The  cyclogenesis  position  was 
to  the  west  of,  and  the  storm  track  over,  an  area  of  maximum 
positive  SST  change.  The  SST  forecast  was  clearly  superior 
to  the  control  run  in  this  case,  especially  considering  that 
ti.e  storm  developed  around  the  fourth  day  of  the  forecast. 

2 .  St  or m  A A  4 

The  tracks  and  SEIS  parameters  for  storm  AA4  are 
si. owl.  in  Figs.  30  and  31.  Storm  AA4  was  forecast  by  both 
the  control  and  SST  runs  at  120  hours  and  was  first  detected 
in  the  analysis  12  hours  later.  The  cyclogenesis  position 
was  in  an  area  of  maximum  positive  SST  change  (1.2  C)  .  3oth 
forecasts  overestimated  the  actual  storm  lifetime  of  72 
hours  by  48  hours.  The  forecast  tracks  were  similar  until 
ti.e  sixth  forecast  (third  day  after  formation),  when  the  SST 
storm  began  erratic  movement.  The  forecast  SEIS  parameters 
in  Fig.  31  are  similar,  with  significant  differences  m 
radius  and  angle  at  the  second  forecast  period  and  radius 
and  ellipticity  at  the  eighth  forecast  period.  The  control 
track  forecast  for  AA4  was  very  sligntly  more  accurate  tnan 
the  SST  track  forecast  during  the  first  three  days.  The 
differences  between  the  forecasts  again  are  mucr.  less  than 
ti.e  departures  of  each  forecast  from  the  analysis. 

'2 .  Storm  AP3 

The  tracks  ana  SEIS  parameters  for  storm  AF3  are 
shown  in  Figs.  32  and  33.  Cyclogenesis  occurred  at  tD  hours 
ir  noth  forecasts  and  ti.e  analysis.  Ti.e  cvcio genesis  ana 


tracks  were  in  an  area  of  positive  SSI  change  (0.4-0.  6  C)  . 
Tne  control  run  forecast  a  weak  storm  witn  a  lifetime  of 
cniy  12  hours  whereas  the  SST  run  predicted  a  lifetime  of  60 
hours.  Since  the  analyzed  storm  lifetime  .as  4d  hours,  the 
SSI  run  more  accurately  predicted  tne  storm  lifetime  and 
track  in  this  case.  The  amplitude  (Fig.  33)  of  this  rela¬ 
tively  weak  storm  was  forecast  rattier  accurately  in  the  SSI 
run . 

C.  SUMMARY 

Significant  differences  in  the  forecast  tracks  ana  SEI3 
parameters  were  noted  in  four  of  the  15  matched  storm 
tracks.  The  control  run  was  more  accurate  ir.  two  of  the 
cases,  correctly  forecasting  November  cyclolysis  wnen  the 
SSI  rur.  did  not  and  more  closely  matching  an  April  storm 

track  wt.er.  the  SST  run  exhibited  erratic  movement.  The  SST 
rut.  was  more  accurate  in  two  April  cases,  forecasting  an 
ana  1 i ze  :  cyclogeaeSxt  wticn  the  control  run  did  not  and 
i-uin turning  a  storm  ..n  wmch  t  ne  control  run  forecast  prema- 

*  lit  cp  clo.ysi...  The  different  method  of  applying  the 

i :  .it  i..  I  oicii.  j  and  tie  increased  vertical  resolution  of 

*  ■  v.  .  s  ;k  IT-  urel  ir.  the  April  runs  made  direct 

*..>»♦. -mi  .  r  » i.  i  A  pci  a  forecasts  difficult. 

.  .  .  ■  .  .  '  w  <  i  •  v  •  r  y  .  r. .-  o  e  r  a  j  1 1.  ;  .  It  appeared  that 

.  . :  i  sit  i  ve  to  tne  effects  of  the 

i ;  ••  i »  [..utiv.ty  translated  into 


FIT.  CONCLUSIONS 


The  introduction  of  time-dependent  sea-surface  tempera¬ 
tures  into  the  NOGAPS  model  resulted  in  significant  changes 
in  mean  and  peak  surface  heat  fluxes  and  precipitation. 
Forecasts  which  use  these  fields  directly,  such  as  model- 
output  statistics,  IGFS  and  NCLAPS-type  models,  may  benefit 
from  the  incorporation  of  time-dependent  sea-surface 
tern pera  tures. 

Differences  between  the  control  and  SST  forecast  storm 
tracks  and  SETS  parameters  were  significant  in  u  of  the  15 
matched  storms.  The  April  SST  run,  which  utilized  the  nine- 
level  H OGA PS  model,  correctly  predicted  cyclogenesis  and 
cyclone  maintenance  for  two  cases  in  which  the  control  run 
did  not.  The  control  run  more  accurately  predicted  a 
November  cvclolysis  event  and  one  April  storm  track.  Only 
minor  differences  in  forecast  positions  and  SEIS  parameters 
were  noted  in  the  other  11  matched  storms.  Differences  in 
the  SEIS  parameters  were  evident  nefore  those  in  storm  posi¬ 
tion.  The  new  nine-level  NOGAPS,  witn  sop histicatic  appli¬ 
cation  of  diabatic  forcing,  seems  sensitive  to  the  effects 
of  time-dependent  SST  ir.  a  way  which  improves  the  accuracy 
of  some  forecasts. 

This  study  was  limited  in  its  application  and  results. 
Further  studies  on  the  response  of  atmospheric  prediction 
*..od  =  ls  to  time-dependent  SST  will  jjc  reguired  to  make 
st at istically  significant  inferences  about  the  effects  of 


suer  boundary  conditions  on  the  guality  of  the  numerical 
forecasts.  Additional  case  studies  witn  the  analyses 
L>:ti  need  to  cover  the  tropical  and  Soutnern  Hemisphere  ocean 
at ca:  srouid  be  conducted.  Additional  model  variables,  suen 
a.  cloudiness,  PEL  thickness  and  static  energy  and  diabatic 


a.r:J  cumulus  heating  should  L>e  analyzed.  Further  attempts 
should  te  made  to  verify  t  ne  heat  flux  and  precipitation 
fe  recasts  of  models  which  incorporate  time-dependent  SSI  to 
determine  model  sensitivity  and  Liases.  The  effect  of  model 
output  changes  on  proposed  long  range  model-output  statis¬ 
tics  and  other  model-derived  forecasts  should  also  be 
eval uated. 


appendix  a 

STORM  TRACK  DATA 

The  11  forecast  storm  tracks  whicn  were  judged  *  o  show 
less  signi  £  icar.  t  differences  are  discussed  it  tuis  appendix. 
I:.  general,  storms  which  were  present  in  the  forecast 
in i t iai  conditions  or  developed  early  in  the  forecast  period 
d^d  not  develop  significant  differences  between  the  control 
and  SSI  forecasts. 

A.  STORMS  DURING  NOVEMBER 

1 .  Storm  N A _1 

The  tracks  and  SEIS  parameters  for  storm  N A 1  are 
mown  in  Figs.  34  and  35.  This  storm  was  present  in  the 
initial  conditions  of  both  forecasts.  Only  minor  differ¬ 
ences  netween  the  SST  and  control  forecast  parameters 
develop  after  six  to  eight  12-hour  forecasts.  Both  forecast 
stores  moved  too  slowly  after  24  hours  and  persisted  longer 
than  the  analyzed  storm. 

2.  Storm  NA2 

The  tracks  and  parameters  for  storm  h'A 2  are  shown  in 
Figs.  3  6-  ar.d  37.  In  this  case  bota  forecasts  accurately 

predicted  storm  lifetime  while  moving  the  storm  too  quickly 
ar.d  ir.  a  more  zonal  path  than  was  analyzed.  Minor  differ¬ 
ences  in  only  two  parameters,  position  and  angle,  can  be 

Sean  at  the  sixth  forecast  period. 

3  .  Storm  N  A  3 

The  traces  and  parameters  for  storm  N A 3  are  shown  in 
rips.  3?  and  3n.  This  storm  formed  at  72  hours  ir.tc  the 


The  tracks  arid  parameters  for  storm  NA4  are  snown  in 
rigs.  4 j  and  41.  This  storm  formed  at  72  hoars.  while  the 
control  and  SST  forecast  parameters  are  similar,  the  control 
storm  lifetime  was  36  hoars  greater  than  that  of  the  SST 
forecast.  The  analysis  and  Loth  forecast  tracks  were 
erratic. 


5  .  Storm  NP1 


The  tracks  and  parameters  for  storm  NP1  are  shown  in 
Figs.  42  and  43.  This  Aleutian  low  was  present  in  the 
initial  condition  of  both  forecasts  and  persisted  for  7  days 
in  each  case.  A  rather  extreme  soutneast ward  deflection  in 
the  analysis  occurred  between  the  third  and  fourth  from  the 
end  positions.  Differences  in  the  SEIS  forecast  parameters 
are  apparent  after  the  seventh  through  ninth  forecasts. 
Only  minor  forecast  storm  position  differences  are  evident. 

f  .  Storm  HP 2 


The  tracks  and  parameters  for  storm  NP2  are  shown  in 
Figs.  44  and  45.  This  storm  formed  at  24  hours  into  the 
ZwD-Lour  forecast  period.  Both  forecasts  underestimate 
storm  lifetime.  He  differences  in  tne  forecast  storm 
tracks,  and  only  minor  differences  in  the  forecast  SEIS 


:  a  r  a  m  et  ers  ,  were  noted. 


S  t  o  r  r  N  P  3 


7  . 

The  tracks  and  parameters  of  storm  NY  3  are  snowi.  i  n 
Fxjs.  46  and  47.  Storm  NP3  was  tna  extra  tropic  a i  transition 
of  a  tropical  cyclone  which  was  not  analyzed  until  43  hours 
into  the  forecast.  SZI5  did  not  trac*.  this  vortex  until  132 
hours  because  the  storm  size  was  less  than  2  grid  lengths, 
loth  forecasts  predict  cyclogenesis  m  the  area  of  the 
actual  extratropicai  transition  approximately  48  hours 
before  xt  actually  occurred. 

B.  ST03HS  DURING  APRIL 

1  .  Storm  A  A_1 

The  tracks  and  parameters  for  storm  AA1  are  snowr.  in 
Frys.  4c  and  49.  This  high  latitude  storm  was  present  in 
the  initial  conditions  of  both  forecasts.  Despite  the  120 
to  132  hour  forecast  storb  lifetime,  only  minor  differences 
in  the  forecast  positions  and  parameters  are  noted.  Tr.e 
analyzed  and  forecast  storm  tracks  are  rougniy  parallel  to 
tne  climatic  ice  edge. 

2  .  Dt_orm  A  A  2 

Ti.e  tracks  and  parameters  for  storm  AA2  are  shown  in 
Figs.  5  0  and  51.  This  storm  formed  12  i.ours  into  the  fore¬ 
cast  period.  Both  forecasts  seriously  overestimated  the 
analyzed  storm  lifetime  by  of  36  hours  i>'j  72  to  34  hours. 

3  .  Btor n  A P_1 

The  tracks  and  parameters  for  storm  A?1  are  shown  in 
I- 1  j  c .  hi  and  53.  This  storm  formed  24  hours  into  the  fore- 
c<;  s  t  period.  Both  forecasts  moved  tne  store,  too  far  i  r.l  an  \ 
an  1  overestimated  ti.e  actual  lifetime  of  2  >  ho  urs  by  4  8 
hours.  Ui(  would  net  exp-.-ct  significant  effects  due  t "> 
tilt  — depei.Uti.t  on  t  n  t  tint  scale  OX  t  n  e  S  e  2  O  r  e  C  a  1'  t  s . 


4.  Ftorn  All 

The  tracks  and  t-arair.eturs  :  jr  store  APT  are  shown  in 
Fiji..  0^  irk  5  5.  This  store,  forte]  at  46  hours  in  the  anal¬ 
ysis  and  toti.  forecasts.  Both  forecasts  under  a  si  ima  te  d  tne 
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periods,  the  a  cpi  i  t  udt  was  too  stall  and  tne  S  SIS  could  no 
loner  track  t  ae  stori. 


Cbtr?  b 


Fiel d- mean  SST  differences  (  C(  during  April  in 
ja)  Pacific  and  (b)  Atlantic  and  during  November 
xn  (c)  Pacific  and  (d)  Atlantic.  Dashed  lines  are 
linear  regression  of  mean  SST  difference. 


Field  statistics  for  sensible  heat  flux  during  November 
for  tiit  Atlantic  region.  (a)  Field-mean  difference 
in  ym-c al/em2-h.  ( c )  Ratio  of  field-mean  difference 

to  field-mean  control.  (c)  Ratio  of  peak  difference 
to  pea)'  control.  (d)  Ratio  of  peak  difference  to 
mean  control. 


4  3 


The  difference  between  the  latent  heat  flux  fields 
of  the  5ST  and  control  runs  at  day  10  for  the  November 
case.  Contour  interval  5  qm-ca l/cm?-h.  Thin  solid 
lines  are  higher  LKF,  thick  solid  line  is  no  change  and 
dashed  lines  are  lower  LHF. 


Field  statistics  for  latent  heat  flux  during  November  in 
in  the  Atlantic  region.  (a)  Field-mean  difference  in 
q c-cal/cn2 - h.  (b)  Eatio  of  field-mean  difference  to 
field-n.ean  control.  (c)  Ratio  of  peak  difference  to  pea 
control.  (d)  Eatio  of  peak  difference  to  mean  control. 


Oi 


Storm  tracks  for  storm  NA5  which  formed  at 
tau  =  72  hoars.  Symbols  mark  12-hour . analysis 
and  forecast  positions.  Solid  line  is  based  on  the 
analyses,  dashed  line  is  for  control  forecast  and 
dotted  line  is  for  SST  forecast. 


{a)  Arj  Lit ude, .  ( l )  angle,  (c)  radius  and 
:or  stern  !i  A  5  rased  on  ar.a  lyses  (solid)  , 
(aashet )  and  SSI  forecast  (Sotted). 


(3)  ellipticit 
control  forecas 


73 


3b.  Similar  to  Fig.  26,  except  for  storm  NA2  which 
was  present  in  the  initial  conditions. 
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